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ABSTRACT 



We investigate the interaction of magnetohydrody- 
namic waves with plasma density inhomogeneities. 
Our numerical study of the full MHD equations 
shows that: (A) Plasma density inhomogeneities are 
a s ource of non-linear generation of transverse com- 



pre isive waves by a plane Alfven wave, and substan- 
tially enhance (by about a factor ot 2) the genera- 
tion of the longitudinal compressive waves. (B) At- 
tained maximal values of the generated transverse 
compressive waves are insensitive to the strength of 
the plasma density inhomogeneity and the initial am- 
plitude of the Alfven wave. (C) Efficiency of the gen- 
eration depends weakly upon the plasma p param- 
eter. The maximum generated amplitude of trans- 
verse compressive wave, up to 32 % of the initial 
Alfven wave amplitude, is reached for about /3 = 0.5. 
The results obtained demonstrate that plasma inho- 
mogeneities enhance the efficiency of the non-linear 
wave coupling. 
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1. INTRODUCTION 



The structure of the solar corona and solar wind is 
governed by the interaction between outward stream- 
ing coronal plasma and the existing strong magnetic 
fields. It is believed that the fast solar wind is as- 
sociated with coronal holes within which the mag- 
netic field is unipolar, open and probably has fur- 
ther fine structuring. The slow wind is believed to 
originate above and at the boundaries of the coro- 
nal streamers where the magnetic fields are bipolar, 
which forces the plasma to flow around regions where 
the field strength is sufficient to slow down the flow. 
The importance of Alfven waves for the energetics 
of the corona an d solar wind has bee n realized since 
their discovery ( Tu fc Marsch, 1995 ). The interest 



in wave phenomena in the solar wind is especially 
timely now, in the context of interpretation of the 
CLUSTER observations. The Ulysses observations 
already revealed a sharp boundary between fast and 



slow wind. Besides, there are indications that there 
is a further fine structuring of plumes in the coro- 
nal holes, so called spaghetti. Therefore, it is of 
prime importance to the study interactions of the 
Alfven waves with these and possibly other struc- 
tures. The weakly non-linear regime, relevant to 
lower cor onal applications, has be e n studied in great 



detail in Nakariakov ct al. (1997 ), Nakariakov et al 
(1998| ), |Botha ct al. (2000|) , |Tsiklauri ct al. (2001|) 



The relevance of this study is prompted by the obser- 
vations of large amplitude Alfven waves in the solar 
wind. In this work we simulate numerically the inter- 
action of a strongly non-linear Alfvenic pulse with a 
one-dimensional, perpendicular to the magnetic field, 
plasma inhomogeneity. 



2. THE MODEL 



We use the equations of ideal MHD 



dV 

p— + p(V -W)V = -Wp 



Of 
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-B x curlB, 



— = curl(U x B), 

divB = 0, 
^+div(pU)=0, 



(1) 

(2) 
(3) 
(4) 



where B is the magnetic field, V is plasma velocity, p 
is plasma mass density, and p is plasma thermal pres- 
sure for which the adiabatic variation law is assumed. 
We solve equations (l)-(4) in Cartesian coordinates 
(x,y,z) and for simplicity assume that there is no 
variation in the y-direction, i.e. (d/dy = 0). The 
equilibrium state is taken to be an inhomogeneous 
plasma of density po(x) and a uniform magnetic field 
Bq in the z-direction. We consider the physical situ- 
ation when initially longitudinal (V z , B z ) and trans- 
verse (V x , B x ) compressive perturbations are absent 
and the initial amplitude of the linearly polarized 
Alfven wave (V y , B y ) is essentially non- linear (rela- 
tive amplitude normalized to the unperturbed value 
A = 0.5). In this case, if the Alfven wave is initially 
a plane wave, the subsequent evolution of the wave, 
due to the difference in local Alfven speed across the 



Figure 1. Snapshot of B y at t = 15.0. Here A = 0.5. 
A = 0.3125, = 2. 



Figure 2. Snapshot ofV x /CA{x) at t — 15.0. Here 
A = 0.5, A = 0.3125, (3 = 2. 



x-coordinate, leads to the distortion of the wave front 
(phase mixing). Hence the appearance of transverse 
(with respect to the applied magnetic field) gradi- 
ents, which grow in time as a power-law. This trans- 
verse gradient is the only source of the transverse 
compressive perturbations, while longitudinal com- 
pressive perturbations are also generated due to non- 
linear coupling to the Alfven wave. 



3. NUMERICAL RESULTS 



As in [Botha et al. (2000| ), fTsiklauri et al. (200l| ), we 
have used the following background density profile 



Po{x) = 3 — 2 tanh(Ax). 



(5) 



Here, A is a free parameter which controls the steep- 
ness of the density profile gradient. In our normal- 



ization , which is the same a s that of Botha et al i 
2000D , [Tsiklauri et al. (200 1[ ), the local Alfven speed' 
is Ca (x) = 1/ \/3 — 2 tanh(Aa;). Full numerical cal- 
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culation of the non-linear set of MHD equations (1)- 
(4), in the above menti oned geometry, has been per- 
formed using Lare2d (Arber et al., 2001). Lare2d 



is a numerical code which operates by taking a La- 
grangian predictor-corrector time step and after each 
Lagrangian step all variables are conservatively re- 
mapped back onto the original Eulerian grid using 
Van Leer gradient limitcr s. This code was also used 



to produce the r esults in Botha et al. (2000), Lsik 
lauri et al. (200l| ). We set up the code in such a way 



that initially longitudinal and transverse compressive 
perturbations as well as density perturbation are ab- 
sent and the initial amplitude of the Alfven pulse 
is sufficiently non-linear, i.e. A=0.5. In the numer- 
ical simulations the Alfven perturbation is initially 
a plane (with respect to ^-coordinate) pulse, which 
has a Gaussian structure in the ^-coordinate, and is 
moving at the local Alfven speed Ca{x): 



B y (x, z, t) = Acxp I — 



{z-C A {x)tf 



(6) 
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Figure 3. Snapshot of V z /Ca{x) at t = 15.0. 
A = 0.5, A = 0.3125, = 2. 



Here 



Here, 5 is a dimensionless free parameter which con- 
trols the width of the initial Gaussian Alfven pulse 
and it was fixed at 0.1 in our numerical runs. Ini- 
tially, V y is non-zero too, and in our normalization 
it is related to B y as V y — —Ca(x) B y , while all 
other quantities are set equal to zero. The simula- 
tion box size is set by the limits —15.0 < x < 15.0 
and —15.0 < z < 15.0 and the pulse starts to move 
from point z — —12.5 towards the positive z's. When 
A-parameter was fixed at 0.125 we had to increase 
our simulation box 1.5 times in order to accommo- 
date fully the density inhomogencity. In this case 
we increased the number of grid point by 1.5 in 
each direction in order to preserve the same sim- 
ulation resolution. Fig. 1 presents a snapshot of 
B y at time t = 15. It can be seen from the graph 
that because of the difference in local Alfven speeds 
(Note, that for x < Ca(x) = l/VE and for x > 
Ca{x) = 1) the initially flat (with respect to x- 
coordinate) Gaussian pulse has been distorted along 
x — axis. This distortion of the pulse front cre- 
ates a transverse gradient (dB y /dx) which is a driv- 



Figure 4- Snapshot of (p — po)/po at t — 10.0. Here 
A = 0.5, A = 0.3125, f3 = 2. 



Figure 5. Dependence of max(\V x (x, z,t)\) versus A 
for /3 = 2 (solid curve) and (3 = 0.5 (dashed curve). 
Here, A = 0.5. 



ing force for the generation of the transverse com- 
pressive wave (V x ). Interesting observation is that 
leftmost (B y (— 15, z, t)) and rightmost (B y (15, z,t)) 
wings of the Alfven pulse decay due to shock forma- 
tion at different rates. Namely, the rightmost wing 
where local Alfven velocity is larger dissipates at a 
larger rate. This can be roughly explained by the 
fact that in the scalar Cohen-Kursrud equation (cf. 
Verwichte et al. (1999 )), which describes the evolu- 
tion of a weakly non-linear Alfven wave until the 
shock formation, the non-linear term is proportional 
to (C\j(C\ — Cg )), where C s is the speed of sound. 
Therefore it is expected that the shock formation 
would be reached quicker where the local Alfven ve- 
locity is greater (in this case rightmost wing). Yet 
another interesting feature can be seen form this 
graph which is appearance of a small bump in the 
solution at x = (phase-mixing region). This could 
be explained by an energy exchange of the Alfven 
wave with the transverse and longitudinal compres- 
sive waves. The latter is expected to happen exactly 
along the x — line where coupling between all three 
modes is strongest. We present solution of Eqs.(l)- 
(4) with the above described equilibrium and the ini- 
tial conditions for the case when plasma-/? is 2. The 
A-parameter was fixed at 0.3125 which corresponds 
to the rather smooth maximum in the dependence 
of max(\V x (x, z,t)\) versus A for this value of (3 (cf. 
Fig. 5 and pertaining discussion below). At time 
t = transverse and longitudinal compressive waves 
as well as density perturbation were absent from the 
system, while at time t = 15 as we see from Figs. 
2 and 3, these physical quantities grew up to a sub- 
stantial fraction of the initial Alfven wave amplitude. 
The transverse compressive wave, V x , is gen erated by 



the B v dB v / dx driving term see for details Tsiklauri 
et ^1. (2001 ). In the homogeneous plasma (A = 0) 
for a plane wave this term is identically zero, thus, 
in the numerical runs with A = we have seen no 
generation of V x at all. This is the plasma inhomo- 
geneity that leads to the generation of the transverse 
compressive wave in this case. The situation with the 
longitudinal compressive wave, V z , is different: even 
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Figure 6. Dependence of max(\V x (x, z,t)\) versus (3. 
Here, A = 0.5 and A = 0.5. 



in the absence of the plasma density inhomogeneity 
the efficient generation of V z still occurs due to the 
BydBy/dz term. By looking at Fig. 3 we gather that 
the solution consists of two parts: first, two positive 
and negative spikes, which dominate the the solu- 
tion by absolute value, which indeed come from the 
phase-mixing region, i.e. part of the solution which 
is caused by density inhomogeneity and the second, 
relatively plane wave front with a substantial deep 
in the middle which is caused by the non-linear cou- 
pling to the Alfven wave. It can be also gathered 
that the second part is moving at about twice of 
the local Alfven speed because (3 = 2. The relative 
density perturbation presented in Fig. 4 has similar 
structure as V z , but it is enhanced in the region of 
interaction. max(\(p — po)/po\) = 0.15 in this case, 
i.e. 30 % of the amplitude of initial Alfven pertur- 
bation. Next, we explore parametric space of the 
problem. In particular, we investigate how the max- 
imal value of the generated transverse compressive 
wave depends on the plasma density inhomogene- 
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4. CONCLUSIONS 
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Figure 7. Dependence of max(\V x (x, z,t)\)/A versus 
A. Here, A = 0.5 and (3 = 1.25. 



ity parameter, A, plasma (3, and initial amplitude 
of the Alfven wave, A. In Fig. 5 we plot depen- 
dence of max(\V x (x, z, t)\) versus A for (3 = 2 (solid 
curve) and (3 = 0.5 (dashed curve). There are two 
noteworthy features in this graph. First, the max- 
imal value of the generated transverse compressive 
wave depends on the plasma density inhomogeneity 
parameter rather weakly, and second, the efficiency 
of the generation of V x is somewhat larger in the 
f3 = 0.5 case than in the case of (3 = 2. In order 
to explore the latter feature further, we made 4 runs 
of Lare2d code for the different values of /3, while A 
was fixed at 0.5. The results are presented in Fig. 6. 
We gather from the graph that there is a given value 
of (3 for which a maximal value of max{\V x {x, z, t)\) 
is observed. Yet another valuable insight can be ob- 
tained by studying the dependence of the maximal 
value of generated transverse compressive wave on 
the initial amplitude of the Alfven wave as our prob- 
lem is essentially non-linear. In Fig. 7 we plot results 
of 4 numerical runs for different values of A, while A 
and [3 where fixed at 0.5 and 1.25 respectively. We 
gather from Fig. 7 that quite unexpectedly the ratio 
max(\V x (x, z, t)\)/A is rather insensitive to the vari- 
ation of the initial amplitude of the Alfven wave, and 
it varies around about 0.27. In other words, irrespec- 
tive of A, the generated transverse compressive wave 
attains about 27 % of the initial amplitude of the 
Alfven wave for this value of (3. Answering the ques- 
tion what determines the observed saturation level 
of maximal value of the generated transverse com- 
pressive wave, we conjecture that, based on the fact 
that the saturation level is rather insensitive to the 
variation of A (see Fig. 5) and it somewhat stronger 
depends of the variation of (3 (see Figs. 5 and 6), 
it is mostly dictated by the given value of (3. Based 
on the information provided by Fig. 6 we conclude 
that the most efficient generation of the transverse 
compressive wave occurs when (3 — 0.5 and, irrespec- 
tive of A and A parameters, transverse compressive 
wave attains about 32 % of the initial amplitude of 
the Alfven wave. 



Our numerical study of the full MHD equations re- 
vealed the following: (A) Plasma density inhomo- 
geneity is a source of non-linear generation of trans- 
verse compressive waves by a plane Alfven wave, 
and it substantially enhances (by about a factor of 
2) the generation of the longitudinal compressive 
waves. (B) Attained maximal values of the generated 
transverse compressive waves are insensitive to the 
strength of the plasma density inhomogeneity and 
the initial amplitude of the Alfven wave. (C) Effi- 
ciency of the generation depends weakly upon the 
plasma /? parameter. The maximum generated am- 
plitude of transverse compressive wave, up to 32 % 
of the initial Alfven wave amplitude, is reached for 
about (3 — 0.5. (D) Inhomogeneity of the medium 
is essential for the generation of weakly compressible 
MHD turbulence. The results obtained demonstrate 
that plasma inhomogeneities enhance the efficiency 
of the non-linear wave coupling. 

It is worthwhile to mention that in this work we 
neglected variation in the y-direction, i.e. we set 
d/dy = 0. This, in turn, prevents linear coupling 
of the Alfven waves to the compressive waves. How- 
ever, validity of this assumption cannot be tested 
observationally so far. Yet another interesting pos- 
sibility of the generation of compressive waves ex- 
ists even in the homogeneous plasma when either 
Alfven wave front is not initially flat with respect 
to x-coordinate or initially flat Alfven wave travels 
over a spatially confined plasma density inhomogene- 
ity. In both cases we should observe effective genera- 
tion of the the compressive waves which possibly do 
not experience suppression caused by the fact that 
d/dx — > oo as in the case considered here and in 



Tsiklauri et al. (200l[ ) 
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